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What is claimed is : 

1. A method for continually or continuously- 
measuring an analyte present in a biological system, 
5 said method comprising: 

(a) transdermally extracting the analyte from the 
biological system using a sampling system that is in 
operative contact with a skin or mucosal surface of said 
biological system; 
10 (b) obtaining a raw signal from the extracted 

analyte, wherein said raw signal is specifically related 
to the analyte; 

(c) performing a calibration step which correlates 
the raw signal obtained in step (b) with a measurement 

15 value indicative of the concentration of analyte present 
in the biological system at the time of extraction; 

(d) repeating steps (a) - (b) to obtain a series of 
measurement values at selected time intervals, wherein 
the sampling system is maintained in operative contact 

20 with the skin or mucosal surface of said biological 

system to provide for a continual or continuous analyte 
measurement ; and 

(e) predicting a measurement value based on the 
series of measurement values using the Mixtures of 

25 Experts algorithm, where the individual experts have a 
linear form 

n 

i=l 

wherein (An) is an analyte of interest, n is the number 
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of experts, An s is the analyte predicted by Expert i; and 

is a parameter, and the individual experts Az^ are 
further defined by the expression shown as Equation (2) 

m 

An, = E a ij P j + z i (2) 

wherein, An; is the analyte predicted by Expert i; P,. is 
one of m parameters, m is typically less than 100; 
are coefficients; and z ± is a constant; and further where 
the weighting value, w if is defined by the formula shown 
as Equation (3) 




where e refers to the exponential function and the d k 
(note that the d ± in the numerator of Equation 3 is one 
of the d k ) are a parameter set analogous to Equation 2 
that is used to determine the weights w ± . The d k are 
given by Equation 4 

where ot jk is a coefficient, P 4 is one of m parameters, 
and where co k is a constant . 

2. The method of claim 1, wherein the analyte is 



97 



0235.04 
2750-0235 .22 
Patent 

extracted from the biological system in step (a) into a 
collection reservoir to obtain a concentration of the 
analyte in said reservoir. 

3. The method of claim 2, wherein the collection 
reservoir is in contact with the skin or mucosal surface 
of the biological system and the analyte is extracted 
using an iontophoretic current applied to said skin or 
mucosal surface. 

4. The method of claim 3, wherein the collection 
reservoir contains an enzyme that reacts with the 
extracted analyte to produce an electrochemically 
detectable signal . 

5. The method of claim 4, wherein the analyte is 
glucose . 

6. The method of claim 5, wherein the enzyme is 
glucose oxidase. 

7. The method of claim 1, wherein the prediction 
of step (e) is carried out using said series of two or 
more measurement values in an algorithm represented by 
the Mixtures of Experts algorithm, where the individual 
experts have a linear form 

BG = w 1 BG 1 + w 2 BG 2 + w 3 BG 3 ( 5 ) 



98 



2750 



0235.04 
-0235.22 
Patent 



wherein (BG) is blood glucose, there are three experts 
(n=3) and BG ± is the analyte predicted by Expert i; w ± is 
a parameter, and the individual Experts BG ± are further 
5 defined by the expression shown as Equations 6, 7, and 8 

BG 1 = p 1 (time) + q x ( active) + r x ( signal) + s^(BG\cp) + t x (6) 
BG 2 = p 2 ( time) + q 2 {active) + r 2 { signal) + s 2 (BG\cp) + t 2 (7) 
SG 3 = p 3 ( time) +g 3 (active) +r 3 (signal) + s 3 (BG\cp) + t 3 (8) 

wherein, BG; is the analyte predicted by Expert i; 
parameters include, time (elapsed time since the 
sampling system was placed in operative contact with 

10 said biological system) , active (active signal) , signal 
(calibrated signal) , and BG/cp (blood glucose value at a 
calibration point) ; p if q if r ± , and s ± are coefficients; 
and t; is a constant; and further where the weighting 
value, w if is defined by the formulas shown as Equations 

15 9, 10, and 11 

w= (9) 

1 d l d 2 d 3 

e 1 + e 2 + e 3 
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d i ^ d 3 



(10) 



d l ^2 d 3 

e + e + e 



(ID 



where e refers to the exponential function and d d is a 
parameter set (analogous to Equations 6, 7, and 8) that 
are used to determine the weights w it given by Equations 
9, 10, and 11, and 
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d x = x 1 {time) + 1^ (active) +y 1 {signal) + 6 1 (BG|qp) + e x (12) 
d 2 = T 2 (time) + P 2 (active) +Y 2 (^i^^^) + 5 2 (BG| cp)- + e 2 (13) 
d 3 = T 3 (time) + (3 3 (active) +Y 3 (signaI) + 5 3 (BG|cp) + e 3 (14) 

where x i7 fjj, Yi and 5 ± are coefficients, and where Si is 
a constant. 

8. The method of claim 1, wherein the prediction 
5 of step (e) is carried out using said series of two or 
more measurement values in an algorithm represented by 
the Mixtures of Experts algorithm, where the individual 
experts have a linear form 

BG=w 1 BG 1 + w 2 BG 2 + w 3 BG 3 (15) 

10 

wherein (BG) is blood glucose, there are three experts 
(n=3) and BG; is the analyte predicted by Expert i; is 
a parameter, and the individual Experts BG ± are further 
defined by the expression shown as Equations 16, 17, and 
15 18 
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BG 1 =p 1 {time c ) + q 1 ( active) + ^(signal) + s 1 (BG\cp) + t x (16) 



BG 2 = Pz ( time c ) + g 2 (active) + r 2 ( signal) + s 2 (BG\cp) + t 2 (17) 



BG 3 -p 3 (tiine c ) + g 3 (active) + r 3 (signal) + s 3 (BG|cp) + t 3 (18) 



wherein, BG, is the analyte predicted by Expert i; 
parameters include, time c (elapsed time from a 
calibration of said sampling system) , active (active 
signal), signal (calibrated signal) , and BG/cp (blood 
glucose value at a calibration point) ; p if q it r if and s t 
are coefficients; and t, is a constant; and further where 
the weighting value, w lt is defined by the formulas shown 
as Equations 19, 20, and 21 



w. = 



1 d, cL d 

+ e 2 + e 3 



(19) 



"i— 2-4-5 (20) 



e 1 + e + e 



e" 3 

„ = e - (21) 

3 d. d, d 

e 1 + e 2 + e 



10 
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where e refers to the exponential function and is a 
parameter set (analogous to Equations 6, 7, and 8) that 
are used to determine the weights w if given by Equations 
19, 20, and 21, and 

d 1 = n 1 {time c ) + $ 1 {act±ve) +y 1 (signal) +b 1 (BG\cp) +e 1 (22) 



d 2 = T 2 (time c ) +P 2 (active) + y 2 (signal) +5 2 (BG\cp) + e 2 (23) 
d 3 = x 3 (time c ) + $ 3 (act±ve) +y 3 {signal) + 5 3 (BG|cp) + e 3 (24) 



where T ix $ if Yi §i are coefficients, and where is 

a constant . 

10 9. The method of either of claim 7 or claim 8, 

which includes further parameters for measurement values 
selected from the group consisting of temperature, 
ionophoretic voltage, and skin conductivity. 

15 10. A method for measuring blood glucose in a 

subject, said method comprising: 

(a) obtaining a raw signal from a sensing 
apparatus, wherein said raw signal is specifically 
related to the glucose detected by the sensing 

2 0 apparatus; 

(b) performing a calibration step which correlates 
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the raw signal obtained in step (a) with a measurement 
value indicative of the subject's blood glucose 
concent rat ion ; 

(c) repeating step (a) to obtain a series of 
5 measurement values at selected time intervals; and 

(d) predicting a measurement value using the 
Mixtures of Experts algorithm, where the individual 
experts have a linear form: 

n 

An = J2 An ± w i (1) 

i=l 

wherein (An) is blood glucose value, n is the number of 
10 experts, An i is the blood glucose value predicted by- 
Expert i; and w r i is a parameter, and the individual 
experts An i are further defined by the expression shown 
as Equation (2 ) 

wherein, Arij, is the blood glucose value predicted by 
15 Expert i; Pj is one of m parameters, it? is typically less 
than 100; a i: j are coefficients; and z i is a constant; and 
further where the weighting value, w it is defined by the 
formula shown as Equation (3) , 



w . - 



(3) 



20 
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where e refers to the exponential function and the d k 
(note that the d A in the numerator of Equation 3 is one 
of the d k ) are a parameter set analogous to Equation 2 
that is used to determine the weights w s . The d k are 
given by Equation 4 



= £ <Vj +(0 * (4) 

where a jk is a coefficient, P j is one of m parameters, 
and where u> k is a constant. 

10 11. The method of claim 10, where in said Mixtures 

of Experts algorithm, the individual experts have a 
linear form 



EG = w 1 BG 1 + w 2 BG 2 + w 3 BG 3 ( 5 ) 



15 wherein (SG) is blood glucose, there are three experts 

(n=3) and BG ± is the analyte predicted by Expert i; w x is 
a parameter, and the individual Experts BG ± are further 
defined by the expression shown as Equations 6, 7, and 8 

BG 1 = p 1 ( time) + q x (active) + r x (signal) + s 1 {BG\cp) + t ± (6) 
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BG 2 =p 2 {tlme) + q 2 ( active) + r 2 ( signal) + s 2 {BG\cp) + t 2 (7) 



BG 3 =p 3 ( time) + g 3 (active) + r 3 (signal) + s 3 (BG\cp) + t 3 (8) 



wherein, EG; is the analyte predicted by Expert i; 
parameters include, time (elapsed time since the 
sampling system was placed in operative contact with 
said biological system) , active (active signal) , signal 
5 (calibrated signal) , and BG/cp (blood glucose value at a 
calibration point) ; p if q if r if and s ± are coefficients; 
and ti is a constant; and further where the weighting 
value, w it is defined by the formulas shown as Equations 
9, 10, and 11 
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+ e 



(9) 




+ e 



(10) 



cL cL d 3 

e + e + e 



(11) 



where e refers to the exponential function and d i is a 
parameter set (analogous to Equations 6, 7, and 8) that 
are used to determine the weights w it given by Equations 
9, 10, and 11, and 

d 1 = x 1 (time) + P 1 (active) +y 1 {signal) +5 1 {BG\cp) + e x (12) 
d 2 = x 2 {time) + P 2 {active) + y 2 {signal) + 5 2 (BG|cp) + e 2 (13) 
d 3 = x 3 (time) + (3 3 (active) + y 3 {signal) + 5 3 (BG|cp) + e 3 (14) 

where x if 3 i7 y^^ and b ± are coefficients, and where is 
a constant . 

12. The method of claim 10, where in said Mixtures 
of Experts algorithm, the individual experts have a 
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linear form 



BG = w 1 BG l + w 2 BG 2 + w 3 BG 3 (15) 



5 wherein (BG) is blood glucose, there are three experts 

(n=3) and BG ± is the analyte predicted by Expert i; w ± is 
a parameter, and the individual Experts BG ± are further 
defined by the expression shown as Equations 16, 17, and 
18 

BG X = p x ( timej + q 1 ( active) + ^(signal) + s l (BG\cp) + t x (16) 
BG 2 =p 2 (time c ) + q 2 {active) + r 2 ( signal) + s 2 (BG\cp) + t 2 (17) 
BG 3 =p 3 ( time c ) + q 3 ( active) + r 3 { signal) + s 3 (BG\cp) + t 3 (18) 

10 wherein, is the analyte predicted by Expert i; 

parameters include, time c (elapsed time from a 
calibration of said sampling system) , active (active 
signal) , signal (calibrated signal) , and BG/cp (blood 
glucose value at a calibration point) ; p it q if r it and s i 

15 are coefficients; and t i is a constant; and further where 
the weighting value, w it is defined by the formulas shown 
as Equations 19, 20, and 21 
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+ e 



(19) 




+ e 



(20) 




+ e 



(21) 



where e refers to the exponential function and di is a 
parameter set (analogous to Equations 6, 7, and 8) that 
are used to determine the weights w if given by Equations 
19, 20, and 21, and 

d 1 = x 1 (tlme c ) +p x (active) + y 1 (signal) +5 x (BG\cp) + e 1 (22) 
d 2 = T 2 (time c ) + p 2 (active) + y 2 (signal) +5 2 (BG\cp) + e 2 (23) 
d 3 = x 3 ( time c ) + 3 3 (active) + y 3 (signal) +5 3 (BG\cp) + e 3 (24) 



where t,, $ if y i and 6 i are coefficients, and where is 
10 a constant. 
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13. ,The method of either claim 11 or claim 12, 
wherein the sensing apparatus is a near-IR spectrometer. 

14. The method of either claim 11 or claim 12, 

5 wherein the sensing means comprises a biosensor having 
an electrochemical sensing element. 



15. A monitoring system for continually or 
continuously measuring an analyte present in a 
10 biological system, said system comprising, in operative 
combination : 

(a) sampling means for continually or continuously 
extracting the analyte from the biological system, 
wherein said sampling means is adapted for extracting 

15 the analyte across a skin or mucosal surface of said 
biological system; 

(b) sensing means in operative contact with the 
analyte extracted by the sampling means, wherein said 
sensing means obtains a raw signal from the extracted 

20 analyte and said raw signal is specifically related to 
the analyte; and 

(c) microprocessor means in operative 
communication with the sampling means and the sensing 
means, wherein said microprocessor means (i) is used to 

2 5 control the sampling means and the sensing means to 
obtain a series of raw signals at selected time 
intervals during a continual or continuous measurement 
period, (ii) correlate the raw signals with measurement 
values indicative of the concentration of analyte 

30 present in the biological system, and (iii) predict a 
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measurement value using the Mixtures of Experts 
algorithm, where the individual experts have a linear 
form 

n 

An = $Z An.w i (1) 

i=l 

wherein (An) is an analyte of interest, n is the number 
5 of experts, Anj^ is the analyte predicted by Expert i; and 
w ± is a parameter, and the individual experts An i are 
further defined by the expression shown as Equation (2) 

m 

wherein, An i is the analyte predicted by Expert i; Pj is 
one of m parameters, m is typically less than 100; a Xj 
10 are coefficients; and z ± is a constant; and further where 
the weighting value, w if is defined by the formula shown 
as Equation (3) 



(3) 



where e refers to the exponential function and the d k 
15 (note that the d £ in the numerator of Equation 3 is one 
of the d k ) are a parameter set analogous to Equation 2 
that is used to determine the weights The d k are 

given by Equation 4 
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where oc jk is a coefficient, P j is one of m parameters, 
and where go* is a constant. 

16. The monitoring system of claim 15, wherein the 
5 sampling means includes one or more collection 

reservoirs for containing the extracted analyte. 

17. The monitoring system of claim 16, wherein the 
sampling means uses an iontophoret ic current to extract 
the analyte from the biological system. 

18. The monitoring system of claim 17, wherein the 
collection reservoir contains an enzyme that reacts with 
the extracted analyte to produce an electrochemically 
detectable signal . 

19. The monitoring system of claim 18, wherein the 
analyte is glucose and the enzyme is glucose oxidase. 

20 20. A monitoring system for measuring blood 

glucose in a subject, said system comprising, in 

operative combination: 

(a) sensing means in operative contact with the 

subject or with a glucose-containing sample extracted 
25 from the subject, wherein said sensing means obtains a 

raw signal specifically related to blood glucose in the 

subject; and 
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(b) microprocessor means in operative 
communication with the sensing means, wherein said 
microprocessor means (i) is used to control the sensing 
means to obtain a series of raw signals at selected time 
intervals, (ii) correlates the raw signals with 
measurement values indicative of the concentration of 
blood glucose present in the subject, and (iii) predicts 
a measurement value at a further time interval using the 
Mixtures of Experts algorithm, where the individual 
experts have a linear form 

BG = w 1 BG 1 + w 2 BG 2 + w 3 BG 3 ( 5 ) 



wherein (BG) is blood glucose, there are three experts 
(n=3) and BG ± is the analyte predicted by Expert i; w ± is 
a parameter, and the individual Experts BG i are further 
defined by the expression shown as Equations 6, 7, and 8 

BG 1 = p 1 {time) +qr 1 (active) + r 1 (signal) + s x (BG\cp) + t ± (6) 



BG 2 = p 2 ( time) + q 2 ( active) + r 2 {signal) + s 2 (BG\cp) + t 2 (7) 



BG 3 = p 3 [time) + g 3 (active) + r 3 ( signal) + s 3 {BG\cp) + t 3 (8) 



wherein, BG ± is the analyte predicted by Expert i; 
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parameters include, time (elapsed time since the 
sampling system was placed in operative contact with 
said biological system) , active (active signal) , signal 
(calibrated signal) , and BG/cp (blood glucose value at a 
calibration point) ; p if q if r if and s ± are coefficients; 
and t ± is a constant; and further where the weighting 
value, w ± , is. defined by the formulas shown as Equations 
9, 10, and 11 



1 ^ d 2 d 3 

e + e + e 



(9) 



(10) 



w 3 = 



d i ^ 



(11) 



where e refers to the exponential function and d ± is a 
parameter set (analogous to Equations 6, 7, and 8) that 
are used to determine the weights w it given by Equations 
9, 10, and 11, and 

d ± = x 1 (time) +p x (active) + y 1 (signal) +5 1 (BG|cp) + e 1 (12) 
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d 2 = i 2 (time) + $ 2 (active) + y 2 {signal) + 5 2 {BG\cp) + e 2 (13) 
d 3 = T 3 (ti/ne) + p 3 (active) + y 3 {signal) +5 3 {BG\cp) + e 3 (14) 

where i ir 3i/ Yi an d 5j are coefficients, and where is 
a constant . 

21. A monitoring system for measuring blood 
5 glucose in a subject, said system comprising, in 
operative combinat ion : 

(a) sensing means in operative contact with the 
subject or with a glucose-containing sample extracted 
from the subject, wherein said sensing means obtains a 

10 raw signal specifically related to blood glucose in the 
subject; and 

(b) microprocessor means in operative 
communication with the sensing means, wherein said 
microprocessor means (i) is used to control the sensing 

15 means to obtain a series of raw signals at selected time 
intervals, (ii) correlates the raw signals with 
measurement values indicative of the concentration of 
blood glucose present in the subject, and (iii) predicts 
a measurement value at a further time interval using the 

20 Mixtures of Experts algorithm, where the individual 
experts have a linear form 
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BG - w 1 BG 1 + w 2 BG 2 + w^BG^ (15) 

wherein (BG) is blood glucose, there are three experts 
(n=3) and BG± is the analyte predicted by Expert i; is 
a parameter, and the individual Experts BG ± are further 
5 defined by the expression shown as Equations 16, 17, and 
18 

BG ± =p ± ( tijne c ) + g x (active) + ^(signal) + s 1 (BG\cp) + t x (16) 
BG 2 = p 2 ( time c ) + g 2 (active) + r 2 (signal) + s 2 (BG|cp) + t 2 (17) 
BG 3 = p 3 ( time c ) + g 3 (active) + r 3 ( signal) + s 3 (BG\cp) + t 3 (18) 

wherein, BG; is the analyte predicted by Expert i; 
parameters include, time c (elapsed time from a 
calibration of said sampling system) , active (active 
10 signal) , signal (calibrated signal) , and BG/cp (blood 

glucose value at a calibration point) ; p if q if r ± , and s ± 
are coefficients; and t x is a constant; and further where 
the weighting value, w it is defined by the formulas shown 
as Equations 19, 20, and 21 

w. = (19) 

1 <*1 <*2 d 3 

15 
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+ e 



(20) 




+ e 



(21) 



where e refers to the exponential function and cL is a 
parameter set (analogous to Equations 6, 7, and 8) that 
are used to determine the weights w if given by Equations 
19, 20, and 21, and 

d x = t x ( time c ) + $ x (actlve) + y x (signal) + 5 1 (BG\cp) + e 1 (22) 
d 2 = x 2 ( time c ) + (3 2 (active) +y 2 (signal) +6 2 (BG|cp) + e 2 (23) 
d 3 = x 3 ( tlme c ) +(3 3 (active) + y 3 (signal) +b 3 (BG\cp) + e 3 (24) 



where x i7 3 i7 and 5 ± are coefficients, and where e £ is 
a constant. 

22. The monitoring system of either claim 20 or 
claim 21, which includes further parameters for raw 
signals selected from the group consisting of 
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temperature, ionophoretic voltage, and skin 
conductivity. 

23. The monitoring system of either claim 20 or 
claim 21, wherein the sensing means comprises a 
biosensor having an electrochemical sensing element. 

24. The monitoring system of either claim 2 0 or 
claim 21, wherein the sensing means comprises a near-IR 
spectrometer . 
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